The approach to scahng in phase-ordering kinetics 
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The influence of the initial fluctuations on the onset of scaling in the quench to zero temperature 
of a two dimensional system with conserved order parameter, is analyzed in detail with and without 
topological defects. We find that the initial fluctuations greatly affect the way scaling is reached, 
while the number of components of the order parameter does not play a significant role. Under 
strong initial fluctuations the preasymptotic linear behavior is replaced by the mean field behavior 
of the large- model, producing the multiscaling to standard-scaling crossover also in the physically 
interesting cases of systems with a small number of components. 
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Much of the interest in phase ordering-kinetics is de- 
voted to scaling behavior of the structure factor . This 
is a late stage phenomenon which takes place after the 
order parameter has reached local equilibrium and the 
typical domain size L(t) dominates all other lengths. The 
structure factor obeys the scaling form 



C{k,t)=L'^{t)F{x) 



(1) 



where x = kL{t), L{t) ~ t~ and, when the order param- 
eter is conserved, z = 3 or z = 4 in the scalar or vectorial 
case, respectively. All what happens before, in the time 
span between the instant of the quench and the begin- 
ning of the late stage, is generically referred to as early 
stage. Usually the early stage theory is identified with 
the Cahn-Hilliard linear theory |^ and early stage studies 
are mainly focused on how and when this theory breaks 
down IQ^. From investigations along this line emerges 
a picture of the early stage as a time regime where quite 
complex and interesting processes take place. In partic- 
ular, it appears that the breakdown of the linear theory 
starts first at short distances and propagates with time 
towards larger length scales Q. 

Recently [H, the early evolution of the phase-ordering 
process with a conserved order parameter has been inves- 
tigated in a different context with the aim of investigating 
the multiscaling to standard-scaling crossover. In the ex- 
act solution of the large- model Q the late stage form 
(0) is replaced by the more general multiscaling behavior 



(2) 



where km{t) is the peak wave vector, x = k/km, (p{x) = 
■>p{x) with ipix) EE 1 - (a;2 - l)^ and kmL(t) - {Int)^/-^. It 
was shown by Bray and Humayu n pQ | (BH) in the frame- 
work of the 1/A^ expansion that (y) holds as a late stage 
property only when A^, the number of the order param- 
eter components, is strictly infinite while the standard- 
scaling form is eventually restored for any finite value 
of A^. The full time evolution of the structure factor for A^ 
large but finite is then expected to display the crossover 
from multiscaling to standard-scaling, with a crossover 



time which grows with A^ and also with the size A of 
the fluctuations in the initial state. Indeed, the exis- 
tence of the crossover together with this dependence of 
the crossover time on A^ and A has been found in the 
numerical solution ||^ of the BH model. 

This result suggests to investigate the early behavior 
of systems with a small number of components and large 
initial fluctuations, for which no early stage theory is 
available. In this Letter we will focus on how the size 
of the fluctuations in the initial condition affects the ap- 
proach to scaling. The influence of the initial fluctuations 
on the structure of the early stage was previously pointed 
out by Elder, Rogers and Desai Q. Here we are able to 
make a detailed analysis of this effect and to complete the 
picture of what happens before scaling is fully developed. 

We consider the zero temperature quenches of a two- 
dimensional system described by the Cahn-Hilliard equa- 
tion 



|^ = -v2[vV-^'(0)] 



(3) 



with = ^_ 1(02)2 r < 0, g > in the two 

cases of Af = 1 (singular topological defects) and A^ = 4 
(absence of topological defects). In the initial state the 
system is uncorrelated 



(x) ■ 0(f' )>= ^<>(.x - x') 



(4) 



and A measures the strength of the fluctuations. We will 
consider the two cases A ^ (small A) and A ~ 



(large A), where 0eg ~ y ~g" equilibrium value 

of the order parameter. We shall flnd that when A is 
large the early stage of the phase-ordering process is ac- 
curately described by the large- A^ model (mean field the- 
ory), both for A^ = 1 and for A^ = 4. Then, the mul- 
tiscaling to standard-scaling crossover is expected to be 
a generic feature of the early stage, when the order pa- 
rameter is conserved and fiuctuations in the initial state 
are large. Furthermore, we shall see that the structure of 
the crossover follows a pattern similar to what one has 
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in the transition from linear to nonlinear regime when 
the initial fluctuations are small . Namely, it turns out 
that also the mean field theory breakdown takes place 
through the propagation of the failure from the short to 
the large length scales. 

The tool which allows to monitor all possible behaviors 
is the measurement of the spectrum of exponents (p{x) in 
over a sequence of time intervals ||]. The form (|^) 
contains ^) as a particular case with Lp{x) = 1. Consid- 
ering C{k,t) in the time interval {Iq^t) beginning at 
and of duration r, after taking the logarithm of Eq.(||) 
in two dimensions we have 



lnC(fc,t) = ip{x)\nL^{t)+\nF{x) 



(5) 



which yields f{x) as the slope of lnC(fc,i) vs InL^(t), 
over the interval {tQ,T), at fixed x. In principle, ^p{x) 
may depend on t^. If this happens scaling does not hold. 
Instead, if ^p{x) does not depend on the time interval 
there is scaling and more specifically there is multiscal- 
ing if </j(x) depends on x, while scaling is standard if ^p{x) 
is independent of x. 

The time evolution of the structure factor has been 
obtained numerically by simulating the equation of mo- 
tion with a computationally efficient cell dynamical sys- 
tem . The details of the computation are the same 
of Ref . 1 1^ . We used values of the parameters such that 
that (plij — 0.9558 and the initial amplitudes of the fluc- 
tuations were A = 2.083 (large A) or A = 8.3 x 10"^ 
(small A). The results were averaged over 30 realiza- 
tions of the initial conditions on a 256 x 256 lattice for 
A'' = 4 (small and large A) and A*" = 1 (small A). For 

= 1 and large A we used a 512 x 512 lattice and av- 
eraged over 90 realizations (up to t — 5432) and over 50 
realizations (up to t = 54177). 

Let us first go over the case of small A. In the be- 
havior of C{k, t) ior N — 1 (Fig. |]) one can recognize 
two sharply separated time regimes, one with k„i con- 
stant and the other with km decreasing with time (see 
also Fig. ^). The latter corresponds to the late stage 
which displays coarsening and a well formed k~^ Porod 
tail for k > km, indicating the existence of sharp in- 
terfaces. For k < km the structure factor displays the 
power law |p^ . Interesting non trivial behavior can 
be observed also in the early stage where km is constant. 
Here two subregimes can be identified. At the very begin- 
ning, as we shall see more precisely below, the evolution 
is linear except for the high wave vector region where an 
incipient Porod tail appears. Later on, but with km still 
time independent, the evolution is nonlinear for all fc's 
due to the appearence of the power law behavior also for 
k < km- The evolution takes place exactly in the same 
way for A'' = 4, except that the Porod tail is replaced 
by exponential decay, due to the absence of topological 
defects g^S). 

Next, in order to probe the scaling behavior, we have 



measured (p{x) in a sequence of seven time intervals 
(Fig. H). The first four curves display a strong depen- 
dence of (p{x) both on x and on the time interval, while 
the remaining three, for t > 500, are essentially flat and 
time independent. Accordingly, the latter ones belong to 
the late stage with standard-scaling already established, 
while the set of the first four clearly illustrates the ab- 
sence of scaling in the early stage. Furthermore, a closer 
analysis of the x dependence of these curves provides ev- 
idence, also in the case with conserved order parameter, 
for the propagation of the breaking of linear evolution 
from high towards low wave vectors, as previously ob- 
served in the nonconserved case In fact, these curves 
display a minimum which moves towards smaller values 
of X as time proceeds and the behavior of (p{x) to the 
left of the minimum can be accounted for on the basis of 
the linearized theory, but not on the right. When A is 
small at the very beginning the non linear term can be 
neglected in Eq.(^) and the analytical form of the struc- 
ture factor can be written as 



C{k,t) = Aexp [2k^ti}{x)] 



(6) 



with km = —r/2 and x = k/km- Then, taking the loga- 
rithm 



lnC(fc,i) 



Ik^t 



V'(x) 



lnL2(t) +lnA 



(7) 



and comparing with Eq.(^ we can identify ^p{x) ~ 

inL™(t) ^(''')- Indeed, the behavior of 1^9(2;) to the left of 
the minimum for curves up to t ^ 200 is well fitted by 
■0(2;) with a growing time dependent prefactor. The be- 
havior of 1^9(2;) for A^ = 4 is identical except for a slight 
change in the vertical scale. 

The picture changes completely when A is large. The 
time evolution of C(fc, t) is displayed in Fig. ^ and the 
main qualitative differences with Fig. ^ are that there is 
no regime with constant km (see also Fig. ^) and that 
it takes much longer to build up the k^ tail for k < km- 
But the real insight into what is going on is obtained from 
the behavior of (p{x). In Fig. ^ and |^ we have plotted the 
time evolution of f{x) over six different time intervals. 
Again ip{x) displays a minimum which moves towards 
lower values of x as time goes on. However, the behav- 
ior is much different from what we had with small A. 
Considering the set of the first four curves for A^ = 1 
(Fig. while on the right of the minimum there is no 
scaling as before, now on the left there is multiscaling 
with ip{x) time independent and well described by ip{x) 
(notice that these curves are magnified by factors of 30 
with respect to those in Fig. For larger times (curves 
from t ^ 1700 onward) the opposite occurs, namely on 
the right of the minimum ip{x) tends to lose the time de- 
pendence as well as the x dependence, while on the left 
there appears a time dependence. This overall behav- 
ior is quite informative. The set of the first four curves 
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shows that in the early stage mean field theory provides 
an adequate treatment of the nonlinearity all the way 
down to a certain length scale (corresponding to the min- 
imum in (p{x)), which grows with time. The description 
of what happens below this cutoff is beyond the reach of 
mean field theory. Later on {t > 1700) mean field the- 
ory breaks down over all length scales, since the appear- 
ance of a flat behavior together with the disappearance 
of the time dependence for x > 1 signals the onset of 
standard-scaling on the shorter distances, while the time 
dependence for a: < 1 shows that the crossover towards 
standard-scaling is under way over larger distances. The 
iV = 4 (Fig. H) system behaves essentially in the same 
way, except that the breakdown of multiscaling for a; < 1 
does not show up in the time of the simulation, indicat- 
ing that the crossover time is larger for = 4 than for 

= 1. It is worth stressing that from the above analysis 
follows that the time regime preceding the fully developed 
standard-scaling under strong initial fluctuations can be 
made to last much longer than the corresponding early 
stage when the initial fluctuations arc small. 

Additional insight into the difference between small A 
and large A can be obtained from the behavior of the 
peak wave vector (Fig. The curves for small A show 
clearly the existence of the early stage with constant 
followed by a sharp transition to power law growth in 
the scaling regime. Instead, for large A the plateau is re- 
placed by a sharp drop of fc~ ^ , soon followed by increasing 
behavior. This initial decreasing behavior of and the 
subsequent change to positive slope are consistent with 
the solution of the large- model Furthermore, the 
asymptotic slope is reached very slowly. We have ana- 
lyzed the growth of k^y^{t) with the power law t^/^ and 
the results for three time intervals are displayed in Table 
1. The interesting feature here is that the asymptotic 
values 1/3 for iV = 1 and 1/4 for iV = 4 are approached 
from above for small A and from below for large A. Thus, 
in particular, for = 1 and A large the effective growth 
exponent before reaching values close to 1 /3 goes through 
a time region where it is closer to 1/4. This is consistent 
with the mean field description of the early stage when 
A is large and shows that the behavior of the effective 
exponent is quite sensitive to the initial condition for a 
long time. 

In conclusion, we have analyzed the approach to scal- 
ing in the quench to zero temperature of a two dimen- 
sional system, which may (A^ = 1) or may not (A^ — 4) 
form stable topological defects, under the effect of weak 
and strong initial fluctuations. We have found that, for 
physically interesting systems, the value of A^ plays no 
significant role, while the strength of the initial fluctu- 
ations is quite important. When A is small there is a 
clearly identifiable and relatively short early stage, with 
no scaling, dominated by the formation of domains on a 
fixed length scale. After these have been formed, coarsen- 
ing begins, the late stage is entered and dynamics obeys 



standard-scaling. When A is large, instead, the initial 
process of domain formation is bypassed and coarsening 
begins right away, but it is not accompanied by standard- 
scaling. Rather, before standard-scaling sets in there is a 
long time transient characterized by multiscaling which 
shows that the early stage evolution is governed by mean 
field theory. So we are now in the position to make an 
assessment of what is the proper place for the large-A^ 
model in the theory of phase-ordering processes for sys- 
tems with a small number of components. It is clear, on 
the basis of the BH result and of more recent work of 
Rojas and Bray that the large- A^ model does not 
give reliable information on the late stage behavior of 
these systems, neither the 1/N expansion offers a viable 
systematic approximation scheme. Instead, as we have 
found, the large-A^ model qualifies for the early stage 
theory, whenever the initial fluctuations are large. 
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FIGURE CAPTIONS 

1. Time evolution of C(fc, t) for small A and A^ = 1. 
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2. Time evolution of (p{x) for small A and = 1. 

3. Time evolution of C{k, t) for large A and 1. 

4. Time evolution of (p{x) for large A and = 1. 

5. Time evolution of (p{x) for large A and = 4. 

6. Time evolution of the peak wave vector. 

TABLE 



Time intervals 


274-1086 


1367-5432 


6837-54177 


A^ = 1 small A 
N = 1 large A 
TV = 4 small A 
TV = 4 large A 


0.47 ±0.03 
0.223±0.008 
0.42 ±0.03 
0.202±0.008 


0.36 ±0.01 
0.275±0.003 
0.295±0.005 
0.220±0.07 


0.351±0.007 
0.322±0.004 
0.266±0.005 
0.235±0.002 



Table 1. Values of the exponent l/z computed by linear 
regression on the data displayed in Fig. ||. 
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